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' We characterize a class of simple FRW models filled by both dark energy and dark matter 

I in notion of a single potential function of the scale factor a(t); t is the cosmological time. 

^ , It is representing potential of fictitious particle - Universe moving in 1-dimensional well 

^y^' V{a) which the positional variable mimics the evolution of the Universe. Then the class of 

^ ' all dark energy models (called a multiverse) can be regarded as a Banach space naturally 

^j^' equipment in the structure of the Sobolev metric. In this paper we explore notion of 

- - metric introduced in the multiverse which measure distance between any two dark 

energy models. If we choose cold dark matter as a reference one then we can find how 
so far apart are different models offering explanation of present accelerating expansion 
phase of the Universe. We consider both models with dark energy (models with the 
generalized Chaplygin gas, models with variable coefficient equation of state wx 



X 
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parameterized by redshift z, models with phantom matter) as well as models basing 
on some modification of the Friedmann equation (Cardassian models, Dvali-Gabadadze- 
Porati brane models). We argue that because observational data still favor the ACDM 
model all reasonable dark energy models should belong to the nearby neighborhood of 
this model. 

Keywords: geometry of dynamical behaviour, dark energy 



1. Introduction dark energy models as a dynamical systems of 
Newtonian type 

Let us assume a class of cosmological models with the Robertson- Walker symmetry 
(i.e. homogeneous and isotropic) of space slices filled by perfect fluid described by 
the general form of the equation of state WcS = = Wcs{a), where Pcs and 
PcS are effective pressure and energy density. Then dynamics of the models under 
consideration, i.e. the FRW models are governed by two basic equations 

^ = > + 3,), (1) 
p = -SH{p + p), (2) 
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where p{a) and p{a) plays a role of the effective energy density and pressure re- 
spectively, a is the scale factor. H = (In a)' is the Hubble function, a dot means 
differentiation with respect to the cosmological time t; SnG = c = 1. 

The first equation is called the acceleration equation or Raychaudhuri equation 
and the second one is the conservation condition. The system of two basic equation 
will be in closed form after postulating the form of equation of state (E.Q.S). We 
assume that Universe is filled by mixture of i noninteracting fluids (in general) such 
that 

Pes = ^Pi = ^ Wt {a) Pi . (3) 

i i 

For example if we assume dust mater (baryonic and dark mater) and dark energy 
X then effective pressure and energy density are function of the scale factor only 



PeS = + wx{a)px{a), 

peS = Pm,oa~^ + px,oa~^ exp 







(4) 



where wxia) = — is a coefficient of E.Q.S for dark energy, 1 + z = a ^ and a = 1 
is corresponding to a present value of the scale factor. 

Hence for assumed form of E.Q.S equation (1) can be written down in the form 
analogous to the Newtonian equation of motion 

In equation 5 V{a) plays the role of the potential function of the particle of unit 

mass moving in the 1-dimensional potential well. Equation (5) reduces problem of 
dynamics cosmological models with dark energy to the problem of classical mechan- 
ics if only energy density satisfies conservation condition (2). 
Of course system (5) admits the first integral in the form 

_ + V{a) = E = const, (6) 

which is usually called the Friedmann first integral. 

It is convenient to rewrite equation (6) to the new (equivalent) form in which 
appears dimensionless density parameters fij = For this aims it is sufficient 
to reparameterize the cosmological time: t T:\HQ\dt = dr. Then we obtain 

yl+v{x) = Uikfl, (7) 

where x = y = a;', ' = ^ (ao is present value of the scale factor), fife = 

-3 

^kfiCi~'^ = , ^ + z = ^ = ^. Without lost degree of generality, k (curva- 
ture index) and density parameter for curvature can be formally incorporated into 
potential function. 
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Therefore the motion of the system is restricted to the distinguished energy level 

E = ^ilfe.o- In the phase {x,y) space trajectories of the system arc determined from 
relation (7). Therefore the classical motion of the system in configuration space is 
restricted to the domain admissible for motion 

Ve = {x:E-V{x) > 0}. 

In the special case of the ACDM model {ilk,o = 0) we have 

Vix) = -^{n^flx-^ + Oa,ox2}, (8) 

V{xiz)) = -i{i^m,0(l + ^) + (1 - ^m,0)(l + ^)"'}, 

where we use constraint relation Qm,o + ^A,o = 1- 

Formally curvature term can be included into potential but then system should 
be considered on zero energy E = level instead of -E = ^^k,o- If we include the 
curvature term into the potential V then an additional additive term appears. 

The phase space {x, y = x') is divided by trajectory of the flat model (f2fe,o = 0) 
y"^ = —2V{x) on domains occupied by close {^kfi < 0) and open models {flk,o > 0). 

In the general case because first integral is valid for any z we obtain constraint 
condition H{z = 0) = Hq in the form 

f^m.O + ^X,0 + ^kfi = 1- (9) 

We complete the different FRW models with dark energy (see Table 1) and different 
generalized FRW models based on modification of the Priedmann equation (see Ta- 
ble 2), appeared in the context of discussion in the origin of the present accelerating 

phase of the Universe [1,2]. For all these models the dynamics is determined from 
the potential function V{x{z)) and corresponding dynamical system is of the form 

x' = y; y' = -^^; where ^ + y(x) EE lf],,o (10) 

i.e. a Newtonian type. Therefore in the phase space (at finite domain) only critical 
points of the saddle or centre types are admissible. It is simple consequence of 
the fact that characteristic equation of the linearization matrix + Vxx{xq) = 
calculated at the critical points {yo = 0, V{xo) = ^^k,o) are admissible. In any case 
they are corresponding to extremum of the potential. If 14 (xq) = and Vxx{xo) < 
then we obtain saddle points (eigenvalues are real of opposite signs). If at the critical 
point Vxx(xq) > then we obtain centre. 

From the physical point of view max is corresponding moment of switching 
decelerate phase dominate by matter on accelerate dominating by dark energy 
effects. 

While potential function V{a) contain all what is needed to characterize the 
system evolution it is possible to reconstruct this modulo Clkfl immediately from 
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case model V{x{z)) relation 
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ACDM model 

non-flat FRW model with A 

FRW model with 2D topo- 
logical defects px = —^Px 

FRW model with phantom 
dark energy px = - ^Px 

FRW model with phantom 
dark energy px = wxPx, 
wx < — 1 fitted 

FRW model with Chaplygin 
gas Px 



FRW model with general- 
ized Chaplygin gas px = 



Px 



A > 0, a = const 



FRW models with dynami- 
cal E.Q.S. parameterized by 
z Px = [wo + wiz)px 

FRW models with dynami- 
cal E.Q.S parameterized by 
scale factor a px = {wo+ 
wi{l -a))px 

FRW model with quantum 
effects origin from massless 
scalar field at low temperar- 
ture (Casimir effect) px = 
— Px,a > 

flat FRW model with dissi- 
pative dark energy PeflF = 
- 3ap"'H- m = -1.5 is 
fixed 



V{x{z)) = -i{n„i,o(l + Z) + {1^ nm,0)(l + Z)-^} 

V{xiz)) 
V{x{z)) 



|{nm,o(i + z) + fifc,o + fitop,o(i + «)"^} 



V{x{z)) = -^{nm,o(l + Z) + nfc,0 + ^phA^ + ^} 

V{x{z)) = -i{n„,,o(l + z) + Qkfi + «ph,o(l + 

V{x{z)) = -\ {n„,,o(l + + + f2ch,o(l + zr\As + (1 - As)(l + 

V[x{z)) = -\ {n^,o(l + + + i^Ch,o(l + z)-\As + (1 - As){\ + ;j)3(l+«)]T^^ } 

V{x{z)) = -\ {n^,o(l + + fifc.o + i^x,o(l + «)3(«'o-«'i+i)-2 exp[3«;i«]} 

V{X{Z)) = -i {n„,,o(l + ^) + fifc.O + f^Jf,0(l + z)3(-o+m+l)-2 e^p[_ 3^]| 
V{x{zy) = -^{nm,o(l -h + fifc.O + fiA,o(l + «)"^ - f^Cass,0(l + zf} 

Vixiz)) = -i {n^,o(l + + ndiss,o(l + z)-^[As + (1 - + z)3(i-'")]T^} 



Table 2. The potential function for 10 cosmological models beyond the standard general relativity 



case model 



V{x{z)) relation 



1 Cardassian type of Friedmann equa- 
tion, £1^ Q = lO"*^ is fixed 

2 Dvali-Gabadadze-Porrati brane mod- 
els (DGP) 

3 Dcffayct-Dvali-Gabadadze brane mod- 
els with A (DDG) 

4 Randall-Sundrum brane models with 

dark radiation and A = 

5 Randall-Sundrum brane models with 
dark radiation and A (RSB) 

6a Shtanov brane models (Branel) 

6b. Shtanov brane models (Brane2) 

7 modified affine gravity (MAG) model 

8 FRW models of nonlinear gravity with 
Lagrangian density proportional to 
Ricci scalar R (NG) 

9 bouncing models with A (BACDM) 

10a models with energy transfer (dark 
matter *-^ vacuum energy sector) 

10b models with energy transfer (dark 
matter *-* phantom dark energy sec- 
tor) 
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the SNIa observations due to the simple relation between the luminosity distance 
and the Hubble function (for flat models) 

-1 



H{z) 



dz \\ + z 



where 



H'^{x{z)) = -2V{x)x- 



(11) 



(12) 



Now wc can define notion of space all dark energy cosmological models in terms of 
space of all dynamical systems of Newtonian type with potential describing dark 
energy model. 

Definition 1. By multiverse of cosmological models with dark energy we understood 
functional space of all dynamical systems of Newtonian type with potential functions 
V{x) ofC^ class. 

Introduced before notion of the multiverse of dark energy models satisfying 
equation of state (3) can be generalized on the case of cosmological models with 
scalar fields. Let us consider for example the FRW model with single minimally 
coupled to gravity scalar field. They can be defined in the form of perfect fiuid 



with energy density 



\(t) + V{(p) and = ^(f) — V{(j)). Therefore potential 



function of the system is U (0) 



iy(0)a^ and then term —j^cp^a^ can 



be shifted into kinetic part because of conventional character of division on kinetic 
and potential parts for the systems with natural Lagrangian. Finally we obtain the 
dynamical system in the form of energy function 



e(a, d, 0, ^) = y - Y^^^a^ + U {a, (/>), 



(13) 



where U{a,(j)) = - ^V{(j))a'^ - \pa'^. 

The analogous construction can be performed for conformally coupled scalar 
field. Therefore if we fix the form of kinetic energy then the potential function 
U (a, 4>) identify quintessence model. 



2. Metric structure multiverse of dark energy cosmological 
systems 

The space of all dynamical systems on the plane can be naturally equipment in 
the structure of Banach space after introducing the Sobolev norm within a space 
of vector field which are tangent to phase curve 

11/11. = max{|/(a;)|,|c^i/|,...,|9./|}, (1) 

where C is closed subspace of the phase space, / is the vector field defined on the 
phase space; / = [x, — fj]"^ in our case. 
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Let US consider metric defined as 5) = ||/— 5 1|. If we consider two models 
of the multiverse then distance between any two dark energy models are 

dco{l,2) = ma^{\{Vi),-{V2U}, (2) 

where the curvature term we shift into potential then system is considered on 
E = level. Therefore two dark energy models are close if their potential function 
and first derivatives are close. The choice of subset C of configuration space can be 
performed in different way. If we are interesting in comparison any two dark energy 
models from view point of present day observational data. At the present epoch 
X = 1 {z = 0) and 1^(1) = — 5 — 5^fe,o and (2) reduces to 



d{l,2) 



1 

2r7= 



max 



C={l}x{y: 0<^n^g»-2y(l)} 



n 



(2) 
m.OJ 



(3) 



+ 



where wc arbitrary assume roughly constraint to ^k,o for example that —0.1 < 
^^fe,o < 0.1, and a set C chosen around flkfl = on the line a; = 1. Of course if we 
introduce the metric as 



del (1,2) 



max 11(^1), 



{V2U, |(Vl)xx-(X^2)x.|} 



(4) 



it will be precisely differentiate between the models. From the observational point 
of view this metric controls both deceleration and jerk parameters. Wc say that 
any two models are close in sense of the metric if both deceleration and jerk 
parameters are close. The C^- metric can be naturally generalized on the case of 
the C metric which can be expressed in dimensionless parameters controlling all 
(r + 1) time derivatives of the scale factor. 

It is convent to rewrite C*^ metric to the new form 



dco{l, 2) = max • 

c 

where we apply relation 



Hi 



(5) 



V^{z = 0) 



dz 



H\z) 



Hi 



1 



and index in the parent denotes kind of model taken from Table 1 or Table 2. 

The form of the metric (5) can be simply generalized on the case of metric, 
namely 



dc^{l,2) = H-^ 



max ■ 
(i,2/)ec 



d_ 

dx 



(Hf,^{x)-Hf,^{x) 



(6) 



^2 (hUx) - Hf,^{x))] _ +4[|(//(^,)(.)-if(%(.)) 



7, 2008 4:8 WSPC/INSTRUCTION FILE art 



8 Marek Szydlowski, Aleksandra Kurek 



For our aims it is more suitable the notion of C'^ norm defined on the class of C*'{E) 
functions, i.e. for / e C''{E), where E is an open subset of M" we define 

(7) 



k = sup + sup | \Df{x)\\ + . . . + sup \\D''f{x) 



The norm 
way 



on the right-hand side of this equation is defined in the following 



\D''f{x)\\ =max 



dxji . . . dxjk 

, the maximum being taken over jl, . . . ,jk = 1, . . . , n. Each of the spaces C''{E) is 
then the Banach space [3]. 

The above metric of relatively character can be useful in comparing to different 
models with the concordance ACDM one. 

Let us illustrate how the Banach metric works for the different dark energy 
models completed in the Table. Let us consider for simplicity the C° metric. Then 
we obtain 



dco{l, 2) = IFi, - V2,\ = n% + (1 - ^'A'fl) 



(1 - n^^h - ^k,o 



Therefore if we take into account independent prior on density parameter for mat- 
ter ri^^o = ^m.o = from independent galactic observation then we obtain 
dco{l, 2) = |rife,o|- The results of calculation of the C° metric for other models (flat 
for simplicity) are 

dco(3,l)=fitop,0 = (l-f^m,o) 

dco(4,l) = f7ph.o = (l-f^^,o) 

dco{5,l)=3\l+Wx\{l-^l*n,,o) 

rfco(6,l) = 3(l-^)(l-n^,o) 

dcoi7,l)=3{l-As){l-n*^,o 
dcoi8a,l)=3{l-n*^^o)\wo + M 
dco{8b,l) =rf(8a,l) 

rfco(9,l)=4(l-OS,_o-^^A,o) 
dco(10,l) = 3(l-^)(l-r^^,o) 

where Wx, -As, wq, Oa.o should be taken from estimation with fixed value of q. 

The metric d describe residuals with respect to the reference ACDM model 
controlled by second derivatives of the scale factor. The higher derivatives of the 
scale factor will be controlled by a more subtle metric: 

dci(l,2) = \{V,Ux = 1) - {V2UX = 1)1 + \{Vi)Ux = 1) - {V2U{x = 1)1 = 



Ha 



+ 



^(iff^W _ +4 [A (iff,, 
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Table 3. Values of ( 



loo 



and dpi foi" models from Table 1 



Case 


dco 


del 


2 


0.00 


0.00 


3 


0.70 


2.10 


4 


0.70 


3.50 


5 


0.60 


2.84 


6 


0.02 


0.08 


7 


0.02 


0.02 


8a 


1.34 


12.16 


8b 


1.22 


13.29 


9 


0.00 


0.00 


10 


0.02 


0.08 



The results of calculation metric for flat models from Table 1 are 





(2,1) = 











dc-i 


(3,1) = 


3(1- 








del 


(4,1) = 


5(1- 


- f^m^o) 






del 


(5,1) = 


3(1- 




3^x1) 


del 


(6,1) = 


3(1- 


-17^,0)(1-^5)(1 + 


3As) 


del 


(7,1) = 


3(1- 


-r!r,,o)(i-^5)(i4 


3Asa) 


del 


(8a, 1) = 


= 3(1 


-r!rn,o)ko + l|{l 


+ 


3ioo(ioo + l)+«;i 

UlQ + l 


del 


(86, 1) . 


= del 


(8a, 1) 






del 


(9,1) = 


8(1- 


- ^m,0 - ^A,o) 






del 


(10,1) = 


= 3(1 


-l]^,o)(l-A5)(l 





In Table 3 we presented values of dco and d(jo for models from Table 1 . The values 
of unknown parameters were estimating from SNIa data (Gold sample [?]) with 
fixed value of q — 0.3. Distances from the ACDM model are also illustrated on 
Figures lA and IB, respectively. 



1.4- 



1,2- 
1,0- 
0,8- 



Fig. 1. A. d^o and B. d^i for flat models from Table 1 
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3. Conclusion 

In the modern observational cosmology appeared the so-called degeneration prob- 
lem manifested by the fact that many cosmological models which predict different 
evolution scenario in the past and in the future becomes still in good agreements 
with observational data. In the classical approach to model selection problem con- 
clusion are made due to the best fit likelihood function for models under considera- 
tion, i.e. so called likelihood ratio test. Unfortunately this method allows to compare 
only nested models therefore it is not appropriate to compare cosmological models. 
In this context the Bayesian framework of model selection may be useful in decision 
which cosmological models is supported by observational data. 

Recent observational data suggest that the coefficient of E.Q.S is very close to 
constant value w ~ — 1 but only w = —1 imply that energy density of dark matter 
(vacuum) is constant. 

In this paper we characterize class of candidates models with dark energy in 
tools of Sobolev metric in the space of all dynamical system on the plane. 

We believe that adequate model with dark energy should lies within the open 
ball at the centre in the ACDM model and with the small radius e. 

Our convince based on the observation that the ACDM model is structurally 
stable, i.e. small changes r.h.s of the system (potential) does not disturb qualitative 
evolutional scenario of the model. Following the Peixoto theorem structurally stable 
dynamical systems on the plane form an open and dense subset in the space of all 
dynamical systems. Therefore structurally stable systems are typical (or generic) 
(structurally unstable becomes exceptional one) and can be approximated by adding 
small perturbations. 

The metric introduced in the Banach space measure how so far we arc from the 
ACDM - concordance model favored by the present observational data. With the 
help of this metric we can compare different proposition in a similar way like with 
the help of the KuUback-Leiblcr metric [5] . 

Our main conclusion is that all simple dark energy models can be successfully 
unified within a scheme of dynamical systems of Newtonian type and potential func- 
tion uniquely characterize different dark energy models. Some of them are generic 
like ACDM model and another exceptional. The models within close neighborhood 
of the ACDM model are typical because its structurally stable. 
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